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ABS TR GW 


A high-resolution, multi-level, primitive equation ocean model is used to examine 
the response of an idealized, flat-bottom, oceanic regime off northern Califorma to 
steady, equatorward, local wind-forcing during the upwelling season. The model has 
open boundaries on all but the eastern coastal boundary on which either free-slip or 
zero-slip boundary conditions are imposed. Time-invariant Winds, either with or 
without a component of wind stress curl, are used as model forcing to spin-up a 
classical two-dimensional, upwelling-induced coastal jet and undercurrent. Since no 
eddies are generated, a stability analysis of the mean flow is conducted which explores 
both the necessary conditions of mixed (barotropic and baroclinic) instability, through 
calculations of potential vorticity, and the sufficient conditions for baroclimic 
instability, through an application of a simple two-layer stability model. Comparisons 
of model results with observations of the coastal jet in the California Current System 
indicate that the location and the horizontal and vertical current shear associated with 
the model coastal jet compare favorably with observations; however the modeled jet ts 
stronger, deeper and wider than the observed jet. Finally, the incluston of wind stress 
curl and the zero-slip boundary condition are demonstrated to be important elements 


in model simulations of the coastal jet. 
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1 INTRODUCTION 


ome OBJECTIVES 

The attainment of sufficient knowledge to adequately simulate observed synoptic- 
miesoscale features of the California Current Svstem (CCS). and to understand the 
Intricate processes responsible for spatial and temporal variations within the system, 1s 
the goal of a numerical modeling effort presently underway at the Naval Postgraduate 
Senool in Vvfonterey, California. The ability to predict positions. intensities, circulation 
patterns, etc. of fronts, undercurrents, eddies and current filaments, given an initial set 
of observations, could provide invaluable assistance to modern naval operations 
Concerned with boundary current acoustic characteristics and associated naval tactics. 

Recently, U.S. Naval battle group tactics have been modified by the predictions 
of a quasigeostrophic (QG) model developed at Harvard University in conjunction with 
the Naval Postgraduate School. The QG model ts initialized by infra-red satellite 
imagery and air dropped bathvthermograph transects across the Gulf Stream and 
prominent Gulf Stream meanders and rings. After initialization, the QG model ts able 
to predict oceanic conditions out to approximately one week. Accuracv is further 
enhanced by dailv model updates and in-situ verification procedures. The initial results 
of this operation look extremely promising and there is potential for expansion to other 
operating areas of the world. With the ever increasing capabilities of high speed 
computers, the possibilitv of using a high resolution Primitive Equation (PE) model in 
m@eeastern boundary current regime, for instance, is feasible. The vorticity signature of 
an eastern boundarv current is generally less than that of a western boundary current 
due to the absence of western intensification, and the observed scale of eddies is also 
smaller. Additionally, while the Gulf Stream, especially its extension region, 1s in part 
ameopen ocean jet and frontal svstem, the CCS is very much associated with the coastal 
upwelling region. This region is characterized by significant horizontal variations in 
Stratification and large (topographically and coastally induced) vertical velocities. For 
@iese reasons, a high resolution PE model mav be more applicable than a QG model 
for use in areas such as the California Current region. Other agencies, as well, could 
greatly benefit from such a predictive model, such as the Coast Guard in search and 
rescue efforts or agencies involved in oil spill recovery, ocean engineering, resource 


exploration, and at-sea waste disposal operations. 


1] 


The intent of this research 1s to provide a contribution towards understanding the 
processes responsible for synoptic-mesoscale variability in the CCS. The tool used in 
this research effort is a ten-level. high resolution PE ocean model with idealized 
topographic capabilitv. The study region extends approximately 300 km offshore from 
the West coast of North America, and it spans the California coastline from Point Sur 
in the south to Cape Blanco in the north (Figure 1.1). The particular realm Gf ieee 
is the effect of wind forcing in the CCS regime. Previous to this effort)an ijpeeem 
baroclinic jet, intending to simulate the mean CCS, was used as the forcing mechanism 
to explore potential instabilities of the mean flow (Batteen er al., 1985). Instability of 
the mean flow occurred in both flat bottom and topographic experiments. It 1s 
believed that wind forcing of the local ocean domain could be an equally important 
process in the CCS and may provide further insight into the dynamics, kinematics, and 


energetics of this complex eastern boundary current. 


B. MOTIVATION AND BACKGROUND 

The question as to how to design a limited series of experiments to study the 
comiplex mechanisms governing the eastern boundary current regime 1s not a trivial 
one. Implementation of wind forcing in a numerical model must be accomplished in a 
logical manner. The exact method and design of wind forcing experiments raises manv 
questions. Scientific foresight and results from past wind forcing studies must be 
intelligently used, especially when confronted with lunited computer resources. 
Procedurally, the most important criteria in the design of wind forcing experiments 1s 
probably to include only one potential generating mechanism at a time, making it 
possible to isolate that effect during the analysis and interpretation phase. 

The elimination of bottom topography in initial wind forcing experiments allows 
the role of wind forcing to be isolated from the possible coupled role of wind forcing 
with bottom topography. The particular developmental stage of the PE model used in 
this study is such that it does not presently include the depth-averaged. or barotropic 
flow, component (although it is in the process of being implemented). It is believed 
that topographic effects would not be well represented unless the depth-averaged 
component is included. This does not. however, dismiss the effect of bottom 
topography as an important mechanism in the CCS. The study of bottom topography 


should be given great consideration in future development of this numerical model. | 
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It is known that local and remote winds drive coastal circulation, the latter by 
exciting coastal waves (Carton, 1984). Wuind-driven coastal upwelling off California 1s 
closely tied to the seasonal cycle because northerly winds along the coast. which are 
favorable for upwelling, are strongest in a climatological sense during summer (Hickev. 
1979). Due to the strength of the local winds and the correlation of local winds with 
upwelling, this study will concentrate solely on the effects of local wind forcing. 

Chelton (1984) points out that the detailed temporal and spatial characteristics 
of the wind field over the California Current are not well known. Probably the most 
complete description of seasonal winds in the California coastal region is given by 
Nelson (1977). who compiled direct ship observations in 1” square areas frompreemnas 
dating back to the mid-nineteenth century. Chelton (1984) additionally notes that the 
spatial structure of the wind field is similar throughout the year but that the magnitude 
changes seasonally. Although the time dependent nature of winds is an important 
consideration, it 1s seen here that an important scale of wind forcing in the CCS region 
is on the order of a season. Numerical studies to date have considered manv different 
time scales of forcing. Carton and Philander (1984) state that “realistic coastal winds 
are an essentially random succession of storms with differing scales. Philandervand 
Yoon (1982) consider periodic winds and Carton (1984) describes the coastal response 
to an isolated storm. This study will focus on the response to forcing at the longest 
time scale considered to be of importance; that being the seasonal time scale. Shorter 
time scales are, indeed, interesting and warrant investigation as well. This study will 
only cover the summer season due to the interest in upwelling mechanisms. 
Furthermore, because the generation period of eddies (less than a month) observed in 
the CCS is short compared to the seasonal time scale, the winds will be input as ume 
invariant (steady). 

Carton and Philander (1984) discuss that at seasonal time scales, two 
characteristics of the wind field have been considered as possibly driving the seasonal 
circulation of the eastern ocean: the alongshore coniponent of wind Stress. anueawie 
stress curl. Bryan and Ripa (1978), Hickey (1979), Chelton (1984), and otherssijaas 
suggested that upwelling due to the curl of the wind stress is an important mechanism 
in establishing pressure gradients. Carton and Philander (1984) conclude that the wind 
Stress curl is important as a forcing function and that it significantly contributes to 
variability near the coast. Anderson and Gill (1975) and Philander and Yoon (1982) 


conclude that the alongshore component of the wind stress. which need not contain 
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any curl, can also contribute to the seasonal variability of the CCS. Chetton (1984) 
points out that the restriction to a zonally uniform wind stress could be an important 
weakness of eastern boundary current models because there are strong cross-shore 
gradients in the wind stress resulting in a nearshore positive wind stress curl. I[t 1s seen 
in the climatological wind record that the winds over the region of the Cahfornia 
Current appear to be strongest approximately 200 km offshore so that there is a 
nearshore positive wind stress curl year round. Additionally, it 1s noted that the 
magnitude of this wind stress curt varies seasonally, approximately in phase with the 
alongshore wind stress. In a recent numerical study on the dynamics of the CCS. 
McCreary er al., (1986) sought to answer the question of what mechanisms force 
equatorward flow off the California coast in a region of positive wind stress curl. An 
important result of that study is that not only the alongshore component of the wind 1s 
responsible for the coastal jet but also that the positive wind stress curl 1s another 
source of equatorward flow. Obviously, wind stress curl plays a crucial role in the 
dynanucs of the CCS and this study will investigate that effect to some extent. 

Another important question which needs to be addressed is whether the sotutions 
are sensitive to the frictional boundary condition, 1.e., free-shp or zero-ship, along the 
rigid eastern boundary of the numerical model. From classical fluid dynamics theory it 
is expected that the water particles in direct contact with a rigid wall will have zero 
tangential velocity. This zero-slip condition has greater implications as the resolution 
Sea Particular model increases. The PE model used in this study has a zonal 
resolution of 8 km and the Zero-shp condition, therefore, is potentially an important 
consideration. This 1s because a uniform tong-shore wind stress having no curl witt 
produce vorticity in the ocean near the wall if the condition of zero-slip is used, but no 
vorticity will be produced if the free-shp condition is used. 

In a study of boundary conditions in homogeneous ocean models, Blandford (1971) 
explored the effect of various boundary conditions (free-slip and zero-slip) in strongly 
nonlinear flow in a barotropic, B-plane model of the ocean. He found that the solution 
changes from a steady flow for free-slip to an unsteady flow with the generation of 
eddies along the western boundary when zero-slip boundary conditions are imposed. 
This study will compare the effects of free-slip versus zero-slip boundary conditions on 


the eastern boundary. 


ie 


With the guidance of important results of past studies, this study of wind forcing 
of the CCS will focus on local. steady wind forcing of a flat bottom regime with the 
additional effect of wind stress curl. Also, both free-slip and zero-slip boundary 
conditions will be explored. The experiments have been designed to proceed in a 
scientifically logical manner with the purpose of investigating one physical mechanism 


aceasulime: 
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feo PE Eawise Kir TION 


Pee BASIC MODEL DESCRIPTION 


1. Model Equations 
The numerical model used in this studv is a ten-level PE model of a baroclinic 
ocean on an {-plane with no depth-averaged flow (the depth-averaged component 1s 
presently being incorporated). The model is based on the hydrostatic and Boussinesy 
approximations with the rigid lid approximation made. Depth of the model is variable, 
when idealized bottom topography 1s used. however. this study will onlv consider a flat 


bottom. The governing equations written in sigma coordinates are as follows: 


a. Momentum Equations: 








du Re co . iczep ob 7 | ee 
ca = aes — : ies ac ra EE 
Py OX Py D Co Ox 
K_ d*u 
— mM + 6.,{u 2.1 
D2 dc dé ) ( ) 
dv -1 dp’ lo dp’ dD 
ie fA 
dt py dy p, D do dy os 
K.. d7s 
=e 5 + 6 4(¥) (2.2) 
Whereo = z,D 
b. Continuity Equation: 
ow = Gu z Cae ; me 
és CX by po 


ee 


c. Vertically Integrated Hydrostatic Equation: 


p = Dj®,Bdé - PLD Baa (2.4) 


ad. Equation of State: 





B= G2 ea (23) 
e. Thermodynamic Equation: 
ee a e ae == = TNA 
dt i Dacc. DGG 
+ O4(t) (2.6) 


In the above equations. all horizontal partial derivatives are on constant sigma 
surfaces. In addition, the variables used are defined in Table 1. Other svmbols in the 
model equations are listed in Table 2, which gives values of constants used throughout 
this study. 
2. Domain Size and Resolution 
The model domain extends 6° in longitude and 6° in latitude (512 km bv 640 
km) and has § km by 10 km horizontal resolution with ten levels in the vertical ime 
is the finest horizontal resolution that is practical for use on the IBM 3033 at the 
Naval Postgraduate School. 
3. Forcing Fields 
The model can either be spun up from rest or initialized with some Kigeam 
current field, e.g., a baroclinic jet representative of the mean Califorma cUrrentaiiee 
addition, a surface heat flux, dependent on the model predicted sea surface 
temperature, can be computed and used to damp the temperature field toward a 
prescribed climatological value. Details are given below. In this study, the model 1s 
spun-up using a representative climatological! wind field. 
4. Fintte Difference Scheme 
In the horizontal. a space-staggered B-scheme (Arakawa and Lamb, 1977) 1s 
used while a sigma coordinate system controls the vertical. The noisefree version of 
the hydrostatic equation in sigma coordinates advocated by Arakawa and Suarez 


(1983) has been implemented. 
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SYMBOL DEIN T iON 





Z height (positive upwards) 

t uume 

p pressure perturbation from a vertical average 

al temiperature 

Lovey eastward.northward. and vertical (sigma) 
velocity components, respectively 

B buoyancy 

bq dynamic adjustment procedure 

5 solar radiation in the ocean 

Woe turbulent vertical heat flux 





5. Boundary Conditions 


The northern, western. and southern boundaries are open using a modified 
version of the radiation boundary condition of Camerlengo and O'Brien (1980). The 
Serer OOUNdaTY, representing the west coast of North America. is closed, and 


Straight, with freeor zero-slip capability. 


Pe SPECIFIC EXPERIMENTAL CONDITIONS 
1. Nlethod of Determining the Wind Forcing 

herclimatolocical study of Nelson (1977) was used to determine the wind 
forcing for the PE model domain, which extends from 36.5N to 42.5N and from 124W 
to 130W (Figure 1.1). Nelson has compiled long-term composite monthly fields of 
surface wind stress on a 1° square area basis along the west coast of the United States. 
Nelson states that the monthly mean data described in the report adequately resolves 
the seasonal cycle, which is the donunant time scale for coastal upwelling (Mooers es 
Geel? 70}). Although the PE model is of a finer horizontal resolution (i.e.; 8 to 10 km) 
than Nelson's analysis (approx. 110 km), it is the best data set available for seasonal 


forcing in the California Current region. 
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TARE 
VALUES OF CONSTANTS Usebep Se. © ere 


NAME 


earth rotation rate 

specific heat of sea water 

drag coefficient 

latent heat of sea water 
constant reference temperamine 
density of air 

density of sea water at T, 
thermal expansion coefficient 
number of levels in vertical 
meridional grid spacing 

zonal grid spacing 

total ocean depth 

latitude of southern boundarv 
latitude of northern boundary 
longitude of eastern boundary 
longitude of western boundary 
time step 

Coriolis parameter 
acceleration of gravity 
biharmionic momentum diffusion 
biharmonic heat diffusion 
vertical eddy viscosity 

vertical eddv conductivity 


surface air pressure 


20 


The choice of selecting a specific span of time over which climatological data 
would be used was driven by interest in the season most favorable for coastal upwelling 
in the California Current System, which has been most often observed to occur during 
the summer (Hickey, 1979). As a result, the months of June, July, and August were 
chosen for the experimental time domain for this study. 

Close inspection of Nelson’s climatological data over the model's region 
during the sunimer months, Figures 2.1 through 2.3, shows a large variability in the 
surface stress field in the offshore direction, and little variabilitv in the alongshore 
direction. As a result. it was decided to use an idealized alongshore wind stress only. 
and to make it a function of the offshore direction only. 

A geographical interpretation problem arises from the fact that the PE model 
has a straight coastline. The cliunatological data. of course, is based on the actual 
position of the California coast, which is not straight. Furthermore, as Figures 2.4 
through 2.6 show, the wind stress curl has a strong dependence on distance from the 
coast. A simple mapping of data from the actual shape of the coast to the straight 
model coast was performed to partly model the effect that the actual coast has on 
values of wind stress curl. 

The reason for exercising care in the preservation of a realistic pattern of wind 
Stress curl is that wind stress curl is the forcing function for the vertically integrated 
mass transport of the wind-driven ocean circulation (Nelson, 1977). fore importantly, 
horizontal variations in mass transport in the surface (Ekman) layer produce surface 
divergence (convergence) corresponding to positive (negative) wind stress curl. The 
Ekman pumping velocity resulting from the divergence is approximately (pf!) times 
the curl of the wind stress (Gill, 1982). Ekman pumping can significantly contribute to, 
and modify, the coastal upwelling regime. 

Another general observation of Nelson’s wind stress data is that the wind 
Stress vectors during the months of June, July, and August, over the model domain, are 
parallel to the coast within about 10-20% accuracy. This means that the magnitude of 
Wind stress vectors represents, to approximately 80-90% accuracy, the alongshore 
Semiponent Of Surface stress. Since the PE model's coastline is oriented due north and 
south, the observed alongshore component is assumed to coincide with the model's 
north component of surface stress. This observation led to the decision to use the 
magnitude of Nelson’s wind stress as the north component of the model stress. The 
zonal (cross-shore) component of the model stress is taken to be zero in all 


experiments. 
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For each month in consideration, a total of 36 data points were extracted 
from Nelson's data set for both wind stress and wind stress curl because the domain of 
the PE model spans 6° of latitude and 6° of longitude and Nelsons data isvoi sau 
resolution. Since the alongshore variation was considered much less significant than 
the cross-shore variation for both wind stress and wind stress curl, an alongshore 
average of the wind stress curl was calculated from Figures 2.4 through 2.6 and then 
averaged over the three months of June, July, and August. The result is shown by the 
circles connected by a dotted line in Figure 2.7. The dominant feature is a region of 
Strong positive curl within 200 km of the eastern boundary, and weaker negative curl 
over the outer 300 km. This time-averaged and alongshore averaged curl reaches its 
maximum Value at the eastern boundary. 

The modeling effort for producing an idealized wind forcing function was 
based on understanding that the wind stress curl 1s a prominent feature which will 
produce a modifying effect on the dynanuics of ocean circulation. The observed three- 
month average, alongshore-averaged wind stress curl, derived fron: Nelson's 
climatology (Figure 2.7), was closely examined for features which could be represented 
by an analytical function. A Gaussian function in the offshore coordinate was found 
to fit the climatological data satisfactorily. The Gaussian function used to represent 


the observed curl ts: 
6z/@x = 0.77exp[-0.5(NR'935.0)7) 11, (2.7) 


where XR is distance offshore in kilometers. Figure 2.7 shows a plot of this analytical 
function (square boxes connected by a solid line) along with the observed values, and 
the fit 1s Seen to be quite good: 

The PE model is designed to use wind input (for convenience Of Windismnes. 
and heat flux coniputations) rather than wind stress curl, but in order to benefit from a 
realistic wind stress curl function and also retain the ability to force the PE Model with 
wind data, idealized wind data were derived from the above Gaussian representation of 
wind stress curl in the following manner. The first step was to calculate the value of 
the wind stress curl (at each u.v gridpoint) using the Gaussian function in Figure 2.7. 
Next, to comipute the surface stress (at the temperature gridpoints) from the wind 
stress curl (at each u,v gridpoint), the wind stress curl was numerically integrated in the 
offshore direction. The constant of integration was obtained from the alongshore and 


three-month average of the data from Nelson (1977). This constant is t= )e0iae 
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dyne/cm? and it represents the value of the alongshore wind stress at the coast. The 
final step in deriving winds from the idealized stress was to utilize the bulk 
aerodynamuc formula to compute winds at model temperature gridpoints. A constant 
drag coefficient of 0.0013 was utilized with an air density of 0.00123 gnicm*. Wind 
values Were defined at model temperature gridpoints because the surface wind speed is 
used in heat budget calculations within the PE model, and therefore it is needed at each 
temperature gridpoint. Lpon completion of heat budget calculations. the PE model 
computes the wind stress and performs a four-way averaging to get the Stress at each 
u.V gridpoint, whereupon it drives model surface currents. 

A plot of the wind stress, obtained by integration of the idealized curl 
function. is shown by the squares in Figure 2.8. It is clear that the proper Shapemams: 
order of magnitude has been captured well. Any difference between modeled and 
climatological wind stress is probably due to the modeling assumption that the curl is 
determined solely bv the offshore variation of the alongshore stress. A plot of the 
alongshore coniponent of the winds computed from the observed and modeled Stress, 
respectively, using the bulk aerodynamic formula, is shown in Figure 2.9. Again, the 
shape and magnitude is well represented and any differences are just reflections of 
differences in the wind stress plots. This idealized wind (having the stress and stress 
curl as shown in Figures 2.8 and 2.7 respectively) is the actual data which was used to 
force the PE model in some of the wind forcing experiments. 

2. PE Model Surface Thermal Damping 

Surface thermal damping becomes a crucial requirement for realistic modeling 
of ocean temperature fields in the presence of atmospheric surface wind forcing. Tlus 
is because the summer heat flux in the California region is a downward [lux duestome 
greatly reduced sensible and latent heat flux produced by the existence of cold 
(upwelled) water. The stronger the alongshore (southward) winds, the greatermage 
upwelling and the colder the surface water and the greater the downward heat flux. 
The downward heat flux is thus negatively correlated with the sea surface temperature: 
hence it is a damping factor. Another way to understand the damping (as contrasted 
to forcing) nature of the surface heat fluxes in the CCS is to recognize that thesia 
downward surface heat fluxes occur where the ocean surface temperatures anemia 
coldest. Heat flux variations related to atmospheric variations in clouds or air 
teniperature are not considered to be an important forcing factor in the CCS. Heat 


fluxes are therefore computed from bulk formulae using prescribed quantities of 
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Figure 2.8 Wind stress. 
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atmospheric solar radiation, clouds, surface air temperature, relative humidity and 
winds (Haney ef al., 1978). The following representative mean values. for the sunimer 


season. are used: fractional cloud coverage of 48°%, solar insolation at top of the 
2 


4 


atmosphere of 983.3 cal'cm*;dav, and surface relative hunudity of 85%. The value of 
solar radiation was obtained from the Smithsonian Meteorological tables (List, 1963), 
the fractional cloud cover was obtained from the satellite cloud atlas of Muller and 
Feddes (1971). and the relative humiditv was computed from climatologica! data of 
Pyvelson (1983). 

Sensible and latent heat fluxes are calculated using exchange coefficients which 
are proportional to the wind-speed-dependent neutral drag coefficient [C4(V)] of Large 
and Pond (1981). Initial conditions for sea surface and air teniperatures, respectively, 
fieme> © Sea SUriace temperature and 15° C air temperature. The sea surface 
temperature used is a representative offshore temperature for the summer, whereas the 
air temperature was chosen so that there would initially be no net heat flux at the 
Suiface with a steady wind speed of ~ 850 cmisec (wind speed used for experiment 1). 
The imposed criteria of no initial heat flux guaranteed that anv heat flux generated 
after model initialization was only in response to the changing sea surface temperature. 
This criteria also assured that the model heat budget acted only to thermally damp the 
eifects of surface wind forcing and not as an additional forcing function. This is not 
only realistic, as noted above, but it also simplifies the interpretation later. 

The damping nature of the computed heat flux is demonstrated best bv 
considering a linearized form of the modeled equations for the total downward surface 
heat flux. Q. which is Q = (ay SSS1Gy TA” is an apparent atmospheric equilibrium 
Bemiperature, SST the mode! surface-layer temperature and y a coupling coefficient 
Strongly dependent on surface wind speed, (Haney, 1985). The corresponding 
relaxation time for the surface temperature is given by Nie Cobh (Game, 21955). 


femere P is the density of sea water, C, the specific heat of sea water, and h the depth 


of the mixed layer. With mean Pe eidcical winds, y~ 25Wm7K7! (Haney, 1971). 
Using h = 50m, a value is obtained for 47! of 100 Gani -Olice —y tS cirectly 
proportional to wind speed, if the wind speed changes on the order of 20°0, the 
temperature relaxation time will also change by 20% vielding a range for 47! of 80 to 
120 davs. This long danyping time will not mask any change in ocean temperature due 


to dynanucal processes which are expected to have a time scale of onlv several davs. 
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Il. RESULTS OF WIND EXPERI 


A, EXPERINIENTS WITH FREE SLIP BOUNDARY CONDI Ties 
1, Experiment | (Uniform Wind) 

In the first experiment conducted. the idealized winds obtained from the 
procedure described in Chapter II], and shown tn Figure 2.9, were a\era@gediimeeae 
cross-shore direction resulting in a constant value (approximately $30 cm sec). The 
forcing for experiment | was, therefore, uniform (alongshore and cross-shore invariant) 
and steady (time invariant). This ideahzed wind was considered to be the mieaes 
complex: vet it still retains some of the realism of observed, long-termij@imie am 
climatological winds during the summer season over the geographical region 
encompassed by the PE model domain. All experiments in this study, unless otherwise 
stated, included surface thermal damping as described in Chapter II. Additionally. the 
eastern boundary (representing the west coast of North America) was set as a free-slip 
condition in this experiment, meaning that the tangential ocean VelocitigS gene 
unrestricted by fricuonal effects. 

At this time, it 1s important to discuss the initial conditions, in particular the 
model temperature stratification for the wind forcing experiments. A pure exponential 
temperature profile, as a function of height (z), was used to initialize all of the 
experiments within this study. This profile has a length-scale of 450 m with a 
temperature at the first sigma level (ocean surface level) of 13° C and a value aipanaas 


depth of 2° C. The exact form of this teinperatureiumenonis asiollg 


T(z) = 2 + 130298. Gar 


The temperature profile was derived from observations used to support the 
Princeton Dvnalvsis model (Blumberg er a/., 1984) and is representative of the long- 
term. mean climatological temperature stratification for the California coastal region. 
Additionally, the Brunt-Vaisala frequency profile, N7(z), was calculated analytically 


using the temperature function in (3.1) from 


\; = Juecine | Gea 


The model temperature profile and resultant N? profile are shown in Figures 3.1 and 
3.2. Since the initial temperature field was horizontally uniform, the initial ocean 
Surments were set to Zero. 

The model was spun up by turning on the wind, at model time zero, to its full 
magnitude. The impulsive wind produced an impulsive surface stress which resulted in 
inertial oscillations of the near-surface ocean currents. These oscillations are 
demonstrated by a time series of the v-component of velocity at a selected gridpoint 
away from the coast (Figure 3.3). The inertial period approximated from Figure 3.3 of 
fev) day compares well with the theoretical inertial period of 2x'f. Off-shore 
Ekman transport was apparent from the negative values of the u-component of 
merocity tliroWshout the model spin-up (see Figure 3.11). Since the wind was northerly. 
and the model domain is in the northern hemisphere, Ekman transport was directed to 
the right of the direction of the surface stress; which is towards the west (represented 
by negative values of the u-component of velocity). The Ekman flow away from the 
rigid eastern boundary produces some very important effects. Water carried seaward, 
in the oceanic boundary layer, produces divergence at the coast. Due to the principle 
of conservation of mass, the void created by coastal divergence must be filled bv 
upwelling of colder water from lower levels. 

(ieseticet o) wpwelline, wihieli brings colder water up from the lower levels. is 
to perturb the initial. horizontally uniform, teniperature structure. The isotherms are 
bent upward near the coast, which establishes a horizontal temperature gradient in the 
upper levels of the water column. The vertical velocities established by upwelling in 
the PE model compare favorably with an analvtical approximation of upwelling rates 
using basic principles and assumptions (see Appendix). Consistent with geostrophy (a 
balance between the pressure gradient force and the Coriolis force) and the resulting 
thermal wind equations {a relation of vertical shear to the horizontal temperature 
gradient), the model spin-up is characterized by a coastally-trapped, alongshore jet in 
the upper 1000 m of the water column. A simple schematic of the coastal jet is shown 
in Figure 3.4 from Gill (1982). The southward flowing jet along the eastern boundarv 
of the model domain will be the focal point of this study. 

Figures 3.5 through 3.7 show contours of the v-component of surface velocity, 
along with contours of temperature. during the spin-up phase at davs 30, 60, and 90 
respectively. The coastal ject develops in response to the increase in the horizontal 


temperature gradient with time. The near-coast, surface, horizontal temperature 
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Figure 3.1 \fodel initial temperature profile. 
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Figure 3.3 V-component inertial oscillations. 
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Figure 3.4 Schematic of coastal jet formation (from Gill, 1982). 


gradient at day 30 (Figure 3.5) is 0.055°/km corresponding to a maximum surface, 
southward alongshore velocity coniponent of 60 cm sec. At day 60 (Figure 3.6), the 
gradient of temperature has increased to 0.063°.km with an associated maximum 
southward velocity of SO cni,sec. Finally, at day 90 (Figure 3.7), the gradveniee 
teniperature has further increased to 0.070°-km with a maxiniuni southward velocity of 
100 cm, sec. 

The vertical structure of the developed jet is best described by coniputing the 
alongshore and ten-day time-average of the v-component of velocity and plotting that 
field as a function of cross-shore direction and depth. Figure 3.8 shows some 
interesting features. The depth of penetration of the southward flowing. coastal jet is 
approximately 1000 m with an offshore extent of approximately 100 Kriya 
miaximium surface velocity hes approxtmately 25 km offshore and the jet axis 1s 
vertically oriented with depth. The first internal Rossby radius of deformation was 
calculated to be 29.3 km. by the method of Feliks (1985). which coincides well witha 
position of the surface maximum, alongshore velocity component in this experiment. 
An undercurrent is present which is characterized as a deep, slow {~ 3 Ciimsaae 
poleward-flowing current underlying the southward flowing jet. This undercurrent has 
a maximum alongshore velocity component at the coast. Additionally, note that the 
weakest vertical shear occurs on the anticyclonic side (offshore or right side) of the 
southward-flowing jet. While it mav not be easy to see from the figure, the largest 
magnitude of relative vorticity also oecurs on the anticyclonic side of the jet in this 
experiment. 

A cross-section of the alongshoreand time-averaged temperature field is shown 
in Figure 3.9. Away from the coast, the initial teniperature stratification of [5° C at 
the surface and 2° C at depth, has only been slightly modified, if at all. However 
coastal stratification demonstrates a classical upwelling regime by the upward bending 
isotherms. As previously described, the horizontal temperature gradients established 
support the geostrophic jet (Figure 3.8). The surface, coastal teniperature has been 
cooled, from an initial temperature of 15° C, to approximately 9° C. This Goan 
teniperature structure has iniportant implications in potential vorticity analvsis and will 
be discussed later. 

2. Experiment 2 (Wind with Curl) 
Experiment 2 differed from experiment 1 only in that the foreing function 


included the important effect of wind stress curl. Idealized wind data (figure sze) 
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derived from the Gaussian representation of wind stress curl (eqn. 2.7), as described in 
Chapter II, was used to force the PE model in experiment 2. The inclusion of wind 
stress curl advanced the complexity, and realism. of model forcing because the wind 
was no longer uniform, as was the case in experiment |. Otherwise, all aspects of the 
first te experiments were identical: 

The idealized wind stress curl (Figure 2.7) has a maximum value at the coast 
of 0.65 dynes’cni- 100 km which decreases to 0.00 dynes'cm*'100 km at a distance of 
187 km offshore. After the zero crossing, the wind stress curl continues touslouam 
decrease to a valuesor -03l1 dvnes/cm~ 100 kni at the western boundary of thems 
model doniain. 

Nelson (1977) described the observed, long-term, mean climatological wind 


stress curl off the California coast as follows: 


Existence of an offshore wind stress maximum results in a line of zero wind stress 
curl approximately parallel to the coast. Positive curl occurs inshore of the 
maximum wind stress. Negative curl in the offshore region is associated with the 
anticyclonic atmospheric circulation over the interior ocean., The positive curl 
near the coast.is related to topography and to local features in the surface wind 
stress distributions. 


The main effect caused by wind stress curl is that the variation of wind stress 
in the offshore direction causes a variation in the Ekman transport which leads to 
convergence and divergence of mass. With the northerly wind blowing along the 
eastern boundary of the model, the positive curl produces oceanic, upper-level. 
horizontal divergence so fluid from below must be “sucked” vertically into the 
boundary laver to replace fluid transported away by mass divergence. The effect is 
Ekman pumping. Negative curl, conversely, creates mass convergence and causes 
downwelling (Figure 3.10). As previously discussed, the wind stress curl goes through 
zero at 187 km offshore and, therefore. the positive wind stress ctirl along the coaam 
causing Ekman pumping, will have a modifving effect on the coastal upwelling process 
described in section II[.A.1. The effect of an offshore varving wind stress curl) asm 
experiment 2, is depicted by contours of the u-component of velocity during the initial 
model spin-up phase. Figure 3.11 clearly demonstrates that Ekman transport is a 
function of the strength of the northerly winds. Referring back to Figure 2.9, it is 
noted that the wind increases offshore to approximately 190 km and Figure 3.11 


displays the corresponding increase in magnitude of the u-component to approximately 
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190 km. The u-comiponent is negative (offshore) everywhere within the domain which 
is also consistent with Ekman transport being directed to the right of the wind. This 
offshore increase in the magnitude of the u-component of velocity is a representation 
of mass divergence resulting in Ekman pumping. Accordingly, proceeding past the 
zero crossing of wind stress curl further offshore, the magnitude of the u-component 
contours begins to decrease which is an indication of mass convergence. 

A cross-section of the alongshore- and time-averaged v-component of velocity 
imemimexperiiment 2 is displayed in Figure 3.12. Striking differences in the structure and 
magnitude of the jet are immediately discernible in this figure when compared with a 
similar plot from experiment | (Figure 5.8). The offshore extent of the surface coastal 
jet In experiment 2 has increased by almost 90 km as compared with experiment 1. 
Note, however, that the effect is very shallow and quite weak. Albeit its relatively 
weak effect, the vertical shear on the offshore (anticyclonic) side of the coastal jet is, 
indeed, strengthened over that in experiment I. The most noticable difference between 
Micwevo proliles is the difference in the strengths of the jet. Figure 3.8 shows the 
Miaimiuni surface velocity of the jet to be -]05 cnisec Whereas Figure 3.12 shows the 
maximum to be only -75 cm;sec. The explanation of this hes with the value of the 
Wind at the coast in each experiment. Experiment | used a constant value of 
approximately -850 cm/sec (domain climatological value) whereas the coastal value of 
Wind in experiment 2 was onlv about -630 cm.sec (coastal climatological average 
fime) Coastal upwelling induced by the wind stress near the coast is clearly the 
dominant factor over Ekman puniping. A difference of onlv 2 m/sec(~ 5 kts) average 
wind speed can obviously produce a dramatic change in the oceanic response adjacent 


to the coast of California when the wind is northerly. 


B. EXPERIMENTS WITH ZERO SLIP BOUNDARY CONDITIONS 
1. Experiment 3 (Uniform Wind) 

The zero-slip boundary condition at the eastern boundary was introduced into 
experiments 3 and 4. The importance of the free-slip versus the zero-slip boundary 
condition was addressed in Chapter [I and it has been the subject of study bv ocean 
circulation modelers since the studies of Bryan (1963) and Veronis (1966). Blandford 
(1971) demonstrated that the large difference between solutions obtained by Veronis 


and Bryan was caused by the difference between the free-slip and zero-slip conditions 
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Figure 3.10 Schematic of Ekman pumping (from Gill, 1977). 
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As in Figure 


Figure 3.12 


used. Veronis’ free-slin condition produced a maxinium velocity at the coast while the 
zero-slip condition used by Bryan produced a maximum velocity offshore. It is Known 
that barotropic instability can occur if the profile curvature changes sign (Haltiner and 
Williams, 1980) and the zero-slip solution of Bryan's certainly satisfied that necessary 
condition. 

Experiment 3, as stated previously, included the zero-slip boundary condition. 
Bwieother aspects of the model were tdentical to experiment |. Figure 3.13 1s a 
normalized plot of the alongshore (v) component of velocity averaged in time (davs 
86-90), in depth (surface to ~ 1000 m). and alongshore. For comparison, a similar 
plot is displayed for experiment | (Figure 3.14). The effects of the zero-shp boundarv 
condition employed tn experiment 3 are immediately apparent. The frictional influence 
of zero-slip on the near-shore tangential velocity (v-coniponent) dramatically reduced 
the coastal velocity as conipared to the coastal velocity generated in experiment | with 
the free-slip condition. Therefore, the main effect of the zero-slip boundary condition 
is to significantly increase the cyclonic horizontal shear on the inshore side of the jet. 
The cyclonic horizontal shear in experiment 1 was .15 x 107! s7! whereas the zero-slip 
condition increased the shear in experiment 3 to .6 xX 107! s7!. Another effect of the 
zero-slip condition in experiment 3 was to move the depth-averaged coastal jet axis 
approximately 8 km farther offshore than that in experiment | (free-slip). | 

As discussed earlier, a major result of Brvan’s zero-slip boundary condition 
study was that the maximum alongshore velocity was located offshore causing the 
necessary condition of barotropic instability (profile curvature changes sign) to be 
satisfied, whereas the free-slip studv by Veronis did not satisfv that necessary condition. 
In this study, the zero-slip condition was not an essential requirement to satisfy the 
necessary condition of barotropic instability. It can be seen that experiment |, which 
had the free-slip condition. also satisfied the necessarv condition (evident bv the jet 
profile curvature changing sign approximately 75 km offshore in Figure 3.14). The 
upwellinginduced coastal jet, simulated in this eastern boundary study, has a velocity 
maximum offshore regardless of the boundary condition used (free-slip or zero-slip). 
In contrast, the western boundary current in the free-slip study by Veronis had a 
maximum velocity at the coast which prevented the profile curvature from changing 
signs. Apparently, the baroclinicity of the modeled upwelling-induced coastal jet 
domunates the positioning of the jet axis over the effects of the frictional eastern 


boundary condition. 
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[tis again interesting to analyze the timeand alongshore-averaged cross- 
section of the v-component of velocity. Figure 3.15 is the profile Genenatcameam 
experiment 3. The effect of a zero-slip boundary condition is immediately apparent in 
this plot. The tangential velocities along the coast in experiment 3 are much less than 
those generated in experiment | (Figure 3.8). Because of this, the coastal isotachs, in 
experiment 3 are sloped upward which significantly reduces the vertical shear on the 
inshore (cyclonic) side of the jet as compared with experiment 1. The result is that the 
jet is more symmetric in the offshore direction. However, the positive Vorlicity omenme 
inshore side of the jet is now about twice as big as the negative VOrlicityeomeeas 
offshore side. This is very different from the pattern in experiment | (Figure sem 
Another interesting result of the zero-slip boundary condition is that the northward 
flowing undercurrent maximum has also moved further offshore along with the coastal 
jec- 

2. Experiment 4 (Wind with Curl) 

As an extension of the zero-slip boundary condition study, experiment 4 
incorporated zero-slip on the eastern boundary along with the wind stress curl as in 
experiment 2. Except for the zero-slip condition, all other aspects of the miodeiiem 
experiment 4 were identical to experiment 2. As would be expected, the basic 
characteristics of the zero-slip boundary condition, as described in the results of 
experiment 3, were reproduced in experiment 4. Again, the vertical shear Ongi@@e 
cyclonic side of the jet was decreased by the upward-bending isotachs near the coast. 
As in experiment 3, this increased the horizontal shear, and hence the relative vorticity 
on the inshore side of the jet. The undercurrent maximum is the narrowest of allyjgmm 
experiments and its alongshore velocity maximum lies furthest offshore duc toma 
combined effects of wind stress curl and zero-slip. Figure 3.16 demonstrates the salient 


features of the mean velocitv profile. 


C. STABILITY ANALYSIS 

It is now well known that, at times, the instantaneous California Current consists 
of “intense meandering current filaments (jets) intermingled with svnoptic-mesoscale 
eddies” (Mlooers and Robinson, 1984). The exact generation mechanism. or source, of 
the observed eddies and their role in the dynamics of the CCS have vet to be 
determuned. One hypothesis currently under investigation is that the eddies are locally 


generated through dynamic instability of the seasonal mean flow. Thomson (1984) 
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16 Same as Figure 


presented a strong case of an observed cyclonic eddy which formed over the 
continental margin of Vancouver Island through the process of barochnic instability 
with an additional. yet secondary, contribution from barotropic instabilitv. He 
furthermore suggested that the energy for the instability was derived from the vertical 
shear of the geostrophic flow associated with a southeastward surface current and a 
northwestward subsurface flow representative of the Cahfornia Undercurrent. 

In the previous sections of this chapter, the coastal jet profile was described for 
each experiment and was shown to satisfy the necessary condition for barotropic 
instabihty. The results of Thomson (1984) suggest that a more important source may 
be barochnic instability of the mean flow. Of the four experiments discussed in this 
chapter, none developed instabihtv and. therefore. eddies and current filaments were 
not generated. However, it is extremely important to analyze the mean flow to 
establish its potential for nuxed (barotropic and barochnic) instability. If potentiai for 
instability indeed exists in the mean flow generated bv the PE model, it mav provide 
insight into the complexities of eddy generation in the CCS. 

Imphcations of the distribution of potential vorticity as a signature of instability 
have been exanuned by Watts and others (Robinson, 1983) for the Gulf Stream. Watts 
stated that potential vorticity, used as a conservative quantity, can serve as a tracer in 
the ocean. Additionally, the cross-stream distribution of potential vorticity may be 
used as a test of predictions by quasigeostrophic theory regarding necessarv conditions 
for unstable growth of fluctuations. The conditions which are necessary for instabihty 
to occur are very useful because, if they are not satisfied, it can be concluded that the 
dynamical constraints will not allow available energy to be released (Gill, 1982). One 
of the conditions is that the cross-stream derivative of potential vorticity must change 
sign somewhere within the domain. Another condition is that the product of the cross- 
Siredin derivative with the basic current must be positive. Also, theory formally 
requires a basic state current which is slowly changing in space and time, so that 
hinearization apphes for perturbation quantities (Robinson, 1983) ; this requirement 1s 
met by the structure of the coastal jets generated in this study. 

The following expression, in cartesian coordinates, used by Watts to examine the 
potential vorticity (q) signature in the Gulf Stream. will by used to exanune potential 


instabihty of the coastal jet generated in this study: 
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Figure 3.17 shows a cross-section of contours of the time-averaged (davs 81-90) 
and meridionally-averaged potential vorticity for experiment | in units of °C m! sec’! 
scaled by 10°. Similarities in the structure of the potential vorticity field between the 
Gulf Stream (Figure 3.18) and the simulated California Current (Figure 3.17) are 
apparent. There is a tendency for potential vorticity to be uniform along isothermal 
surfaces and also to change vertically in the same manner as the offshore temperature 
stratification. The potential vorticity field for experiment 1 1s further characterized by 
the following features: 

e The range of potential vorticity lies between 0.0-2.5 X 10° °@ im! sec™ 


e <A relative nuninmwum exists in the surface laver caused bv weak stratification due 
to turbulent Vertical nuxing. 


e <A stronger relative nuninium exists in a thicker surface laver near-shore due to 
upWelling induced weak stratification. 


e <A relative maximum in se seianeh vorticitv exists in the middle to western part of 
the domain at a depth of ~ 90 m (the “Seasonal” thermocline in the model). 


e A change in the sign of the slope of potential vorticity contours is apparent in 
the east-west direction of the upper levels, near-shore. 


One difference between the potential vorticity field in the Gulf Stream and that in 
the simulated CCS is the difference in magnitude of potential vorticity. The magnitude 
of potential vorticity in the Gulf Stream is approximately twice that in the simulated 
CCS. This difference is mainly attributable to the difference in temperature 
stratification between the Gulf Stream: and the CCS. From the isotherms depicted in 
Figure 3.18 it can be seen that the strongest stratification in the “Gulf “Sipe 
approaches ~ 10° C over 140 m or 0.07 °C, m. whereas from Figure 3.9 the maximum 
temperature stratification for experiment | is only ~ 3.5 “C over 140 m or 0:03) Gia 
Note that Figure 5.18 was constructed by the approximation of potential vorticity, q 
~ f dT'dz, whereas Figure 3.17 was constructed from the full form of potential 
vorticity (eqn. 3.5). The neglected term in Watts’ approximation was found to 
contribute between 10 and 20°% of the potential vorticity in most Cases saci 


particular, the second term donunated in experiment | on the near-shore, cyclonic side 
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of the coastal jet. The scaled form of potential vorticity was used in the Gulf Stream 
analysis due to limitations in data. As a result of this, only qualitative comparisons 
fave been made between analvses. 

As previously discussed, a necessarv condition for instabilitv is that the cross- 
stream derivative of potential vorticity must change sign somewhere within the domain. 
To examine the potential vorticity field for this condition, the horizontal gradient was 
computed (€q 0x) and multiplied by one grid length (AN) for scaling purposes. The 
resultant cross-stream derivative is displaved in Figure 3.19 for experiment 1. As 
expected from Figure 3.17, the cross-stream derivative changes sign in the upper levels 
across the extent of the simulated coastal jet. The necessary condition for instability of 
the mean flow was therefore, easily satisfied in experiment I. 

The other experiments conducted in this study displaved very similar patterns of 
the cross-stream gradient of potential vorticity and, therefore, all will not be shown. 
All of the experiments satisfied the necessary condition for instability. There was no 
significant difference in the gradient of potential vorticity between the experiments with 
Uniform wind and those which included wind stress curl. However. the gradient of 
potential vorticity in the experiments with Zero-slip was stronger (both positive and 
negative) than those with the free-slip boundary conditions because of the additional 
Somiioulion from horizontal shear as discussed in section IJ1.B.1. (see Figure 3.20). 

Why. then, if the necessary condition for instabilitv has been satisfied in each 
experiment, has the mean flow not gone unstable? The answer to that question is not 
readily apparent. [t is useful, though, to explore a simple barochinic instability model 
to deternune what the suflicient conditions are for baroclinic instability of the mean 
flow. 

The simplest model which can incorporate baroclinic processes is a two-layer 
foogde! Following the method of Holton (1979), the ocean is divided into two discrete 
lavers. In all of the experiments conducted in this study, the southward-flowing coastal 
jet extended to approximately 1000 m in the vertical. For this reason, the upper model 
layer extends from the surface to 500 m depth and the lower layer extends from 500 to 
F000 m. See Table 3 for the arrangement of variables in the vertical for the two-laver 
baroclinic model. 

The QG vorticity equation is applied at levels 1 and 3 whereas the QG 
thermodynamuc equation is applied at level 2. The streamfunction, y, is obtained by 
linearly interpolating between -250 m and -750 m. This forms a closed set of prediction 


equations with B = 0, as in the PE model. 
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The length scale used in this model is the Rossby radius of deformation (47!) 
defined bv: 


7 
f 


Ne —— , 
C( A778 


(3.5) 


Where Az is the layer thickness (500 m) and Go is the stability parameter defined by: 
o=agoT dz. (3:6) 


The phase speed (c) is obtained from solving the closed set of QG equations, by 


assumung Wave-type solutions, yielding: 





kK? = 2h Ja 
c= YM = OY eae ae (29 
where the vertically averaged current (V,y,) is defined as: 
Via ; 
V Vf = . ; (3.8) 


64 


and the basic state “thermal wind” (V7) is defined as: 


Vr oe (3.9) 
; see : ; 4 se : 
For waves with meridional wavenumbers (k) satisfying k~ < 2A~, there exists an 
imaginary component of the phase speed (c.). This criteria defines the critical 


wavelength (L.) as follows: 


L. = 22 (mr) . (3.10) 


Cc 
OT 


I a ee || ae (3.11) 


Since B = 0, the criteria for instability does not depend on the magnitude of the basic 
state “thermal wind” (V7). All wavelengths longer than L, are unstable even for very 
small vertical shear. 

Any wave longer than this critical wavelength will amplifv at the growth rate (@) 


where 


The e-folding growth time is simply the inverse of @. Note that @ depends on V+ and 
feed is a fUnction of the meridional wavenumber, k. It is a maximum at an 
intermediate wavenumber (0 < k*< We 

Table 4 shows the inverse Rossby radius of deformation, the critical wavelength 
(L.), the basic state “thermal wind”, and the most unstable wavelength (L_.) with its 
associated e-folding time for experiments 1-4 and for an imposed, unstable baroclinic 
femexperiment (BCJET) from a previous flat bottom study (Batteen ef al., 1985). In 
this analysis, the computations were made using the time (days 81-90), zonal (within 
150 km of the eastern boundary), and meridionally averaged data from the PE model 
experiments. 

The e-folding times for all of the experiments, at the most unstable wavelengths, 
are very comparable. The computed values of the e-folding times are only an 
approximation to the true baroclinic growth rate due to the simple two-layer baroclinic 
instability model used, and therefore the difference in e-folding time between all of the 


experiments 1s probably negligible. 














TABLE 4 
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The short e-folding time, however, would indicate that eddy development should 
be apparent after about a week or so of model integration .once this unstable state had 
been achieved. Indeed. eddy generation was apparent in the experiment with the 
imposed baroclinic jet after about 25 days of model integration. Notice that there is a 
large difference in the critical wavelength and associated Rossby radius of deformation 
between experiments 1-4 and the imposed baroclinic jet experiment. Possibly the 
reason eddies have not been generated in the wind-forced experiments lies in the 
difference between critical wavelengths. 

Figure 3.21 is a plot of the e-folding times as a function of meridional wavelength 
(L = 2n k) for the imposed baroclinic jet experiment and experiments 1-4. The two- 
laver baroclinic model was applied over a ten-day averaged period of model days $1-90. 
The difference in critical wavelengths is readily apparent in the figure. Additionally, 
the most unstable wavelength is ~ 80 km in the imposed baroclinic jet experiment 
versus ~ 115 km in the windforced experiments. 

To deternune the baroclinicity of the experiments throughout the spin-up period, 
ten-day averages were computed and the two-level baroclinic model was applied for 
three different periods in experiment |. Figure 3.22? shows this progression andgyimem 
clear that the baroclinic flow in this experiment is becoming more unstable throughout 
model spin-up. At days 21-30 the shorter wavelengths had a definite preferred growth 


rate while at davs $1-90, the longer wavelengths are not considerably less unstable than 
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the shorter wavelengths. Supemmposed on Figure 3.22 is the familiar result from the 
imposed baroclinic jet experiment. This figure seemis to imply that if the exoerimicame 
were continued eddies may be generated. The problem with model continuation is that 
the desired realism would be lost. Re-design of model spin-up would be required 
including time variant wind forcing and variable parameterization of surface heat fluxes 
to simulate conditions indicative of the fall season. Possibly, the best solution would 
be to simulate an annual cycle to spin the model up. Additional computer resources 
would be required for an experiment of that scale. 

Once a stability analvsis of the mean flow is completed, it is necessary to examine 
model heat and momentum diffusion and its associated damping (e-folding) time. If 
the damping time is on the order of the growth rate for the unstable wavelengths, 
instability of the mean flow will be suppressed. 

The PE model used in this study incorporates biharmonic heat and momentum 
diffusion. Following Holland and Batteen (1986), the QG baroclinic mode vorticity 


and thermal equation has the form: 
(Vo -k-) YU, = Avo y > he (3.13) 


Where w is QG “temperature” (YW ,-W, in two-layer model), A is biharmonic eddy 
viscosity, B is biharmonic eddy diffusion, A is inverse Rossby radius (as before), and 
only the damping terms are kept in (3.13). 


Assunung wave numbers (k, lL), (3.13) becomes: 


A(k? + 2) + 4°B(k? + 1°) 
ko + 1 + WP 


or w=ywoel . (3.15) 
The damping (e-folding) time of baroclinic modes is therefore 


; (ke + FP + 22) / 
v = = 53 5 > as : (sale) 
Ak? = 17)? bo ele 


Table 5 shows damping rates (7) and time scales y for various wavelengths 


‘las in the PE model. It can be seen that 


using values of A = B = 2 < 10)” ences 
for the most unstable wavelengths predicted by the two-layer baroclinic model for 


experiments 1-4 (~ 100 km), the damping time scale is ~ 95 days: SinGouiine 
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baroclinic instability e-folding time was previously shown to be ~ 3-4 days, diffusive 
damping would appear to be negligible at that wavelength compared with the 
baroclinic growth rate. For longer wavelengths, which also have small e-folding umes 


due to baroclinic instability, the diffusive damping ts entirely negligible. 
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PeeovirAnIsON OF MODEL RESULIS Will OBSERVATIONS 


A. HISTORICAL NOTES ON THE CALIFORNIA CURRENT SYSTEM 

A review and synthesis of existing literature on the CCS through 1978 1s 
presented by Hickey (1979). Throughout the history of observations in the CCS, a set 
of names have been assigned to the various poleward and equatorward flows of the 
CCS. The California Current is the equatorward flow; the California Undercurrent 
refers to poleward flow over the slope; the Davidson Current is the poleward flow 
north of Point Conception during fall and winter: and the Southern California 
Countercurrent refers to the poleward flow south of Point Conception in the California 
Bight (Hickey, 1979). Hickey (1979) points out that as our sampling resolution 
increases, confusion in terminology arises and that: “The confusion in nomenclature 
parallels a growing uncertainty about the temporal and spatial variability of the various 
currents and the dynamical relationships among them.” The decade following 1978 has 
Vet to provide clear-cut definitions of the structure of the CCS and has only helped to 
complicate the picture by demonstrating that the CCS 1s more variable, both spatially 
and temporally, than previously thought. 

Regardless of the complex variabilitv, the flow which has been denoted the 
California Current generally has both an offshore and a nearshore maximum in its 
alongshore component. The large scale equatorward flow is believed to be fed by the 
northern limb of the North Pacific gvre (the West Wind Drift). This large scale 
alongshore component of southward flow has been observed to exceed distances of 500 
km offshore in some locations. The nearshore maximum, or coastal jet, escaped 
attention in early studies of the CCS primarily because the majority of hydrographic 
data had station spacing of about 50 km. Recently, however, the emphasis has shifted 
to conducting surveys of finer resolution in an attempt to resolve the synoptic- 
mesoscale structure of the coastal jet and its associated meanders, seaward jets, and 
eddies. 

In this study of wind forcing experiunents in the CCS, the large-scale component 
of equatorward flow associated with the North Pacific Gyre has not been included. 
The experiments were conducted with an open boundary regional model with local 


wind forcing and for clarity of interpretation, large scale gyre circulation was neglected. 
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In future developments of this model, it may become necessary to include large scale 
gvre circulation in the quest for realistic simulations of the CCS. For the above 
reasons, Only observations of the coastal jet and its associated undercurrent will be 


used for comparison with model results. 


B. THE COASTAL JET 

In Chapter III, the model spin-up phase was described in detain 
description was based on classical, two-dimensional coastal upwelling processes. The 
resultant coastal jet has been a prominent feature in descriptions of shelf circulation in 
the Pacific Northwest during the upwelling season (Hickey. 1979). Actual descriptions 
of the coastal jet, however. indicates that it varies widelv in space and time. Recently, 
hydrographic and XBT surveys have demonstrated that instantaneous cross-sections of 
the California coastal region show great variability. Kosro (1986) concluded that 
synoptic maps of near-surface currents off northern California deviated substantially 
from the classical description of two-dimensional wind driven upwelling. Huyer and 
Kosro (1986) conducted six survevs near Point Arena, California during the upwelling 
season in 1981 and 1982. The simple two-dimensional. equatorward baroclinic coastal 
jet was observed during only one of the six survevs (20-22 April 1982). The simple 
two-dimensional coastal jet is almost invariably resolved in the averaged velocity fields 
off California during the upwelling season. Huver and Kosro (1986) found that the 
average velocity fields (computed over all 6 surveys) consistently showed baroclinic, 
equatorward flow in the upper layers and weak poleward flow along the shelf-break at 
a depth of ~ 150 m. Additionally, it was noted that the eddies and meanders which 
dominated the velocity field during most surveys were absent from the average field. 
All of their surveys revealed a baroclinic equatorward jet which wound around eddies 
and at times turned directly offshore. It is hypothesized that since the classical form of 
the coastal jet appears in the average velocity fields. but only rarely in @ Svmoime 
survev, the upwelling induced coastal jet becomes unstable off northern California 
during the upwelling season. 

The classical two-dimensional coastal jet 1s somewhat more prevalent in its pure 
form off the coast of Oregon. In two successive vears (1965 and 1966), the igi 
zone of the coastal upwelling region off Oregon was studied by Mooers ef al. (1976). 
They found that the coniputed velocity field showed equatorward surface flow and a 


poleward undercurrent below the inclined frontal layer. Both the southward surface 


—~] 
ty 


flow and the northward undercurrent exhibited a “jet-like structure.” The southward 
jet was most intense (~ 20 cms) within 20 km from the shore and in the upper 40 m. 
The poleward undercurrent was most intense (~ 10 cmi’s) between [0 and 40 km 
offshore and at a depth of $0 to 200 m. 

Huyer and Kosro (1986) state that, as previously discussed, a baroclinic coastal 
jet was manifested in classical form during their 20-22 April 1982 survey off northern 
@amiormia. ihe strength and position of the jets core (over the shelf-break) were 
nearly uniform over the entire alongshore extent of the survey region (~ 100 km). 
This anisotropic characteristic (miore rapid variation cross-shore than alongshore) 1s a 
feemmrent feature Of the coastal jet im iis Stable form. The coastal jet was further 
Giaracterized bv a clearly defined core located about 25 km from the coast with 
Peminward velocities = 50 cms at 20 mi depth. Increasing in depth to 80 m, the core 
was located farther from shore (~ 40 km) with maximum southward velocities less 
than 30 ni’s. Available data suggested that there existed a narrow (10-20 km) poleward 
Mudercurrent at ~ 150 m depth with maximum velocities of about 15 cm/s. 

Strong horizontal and vertical shear is a persistent characteristic of the baroclinic 
coastal jet. Huver (1983) found the maximum Vertical shear to be relatively constant 
(~ 4 x 10° s!) during the upwelling season in spite of significant wind stress 
fluctuations. Huyer and Kosro (1986) used current meter data to deternune that 


average maximiuni Vertical shear was 4-5 X Ore 


Mm tdeupper so Mm) over tic-outer 
shelf during the 1981 and 1982 surveys. Kosro (1986) computed vertical shear to be 
greater than 2 X 10°? s*! in the coastal jet off northern California during CODE with 
observed cyclonic horizontal shear exceeding 10°* s"! which is the local value of the 
Coriolis parameter (f). Mooers er al. (1976) computed mean geostrophic vertical shear 
from the thermal wind equation to be ~ 4 x 107 s7! off Oregon which compared well 
with their measured current shear of 4.6 xX 107° s in the coastal [cmon Orecon. 
Kosro and Huver (1986) collected data across a seaward tet off northern California and 
suggested that it was continuous with the coastal jet. Horizontal shear on the cyclonic 
side of the jet was estimated to be on the order of 107 s! (f) and on the order of f'3 or 


fess (Or the anticyclonic side of the jet. 


= COMPARISONS 
The most important consideration in a comparison of model results with 


observations is which model assumptions might cause the model results to deviate most 


8: 


sigmificantly from nature. In this study there are two rather artificial assumptions, the 
flat-bottom and the straight coastline, which could significantly alter the results from 
that observed in nature. As discussed in Chapter I], however, the flat-bottom and 
straight coastline are crucial in interpretation of early modeling results. It is with this 
in mind that the following comparisons are presented. Table 6, below, is a synopsis of 
some of the more interesting features of the coastal jet generated in experiments 1-4 
along with generalized features of a classical. upwelling-induced, two-dimensional 


coastal jet occasionally found off the coast of northern California and Oregon. 
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The letter designations in Table 6 have the following meanings: 
A. Maxunum jet velocity. 

Offshore location of surface jet axis. 

Offshore location of depth averaged jet axis. 

Oiisiene extent of jer. 

Deprhros jet. 

Maximum undercurrent velocity. 

Width of undercurrent. 

Maximum vertical shear (x 107°). 

Maximum horizontal cyclonic shear. 


Maximum horizontal anticvclonic shear. 


1 


Za 2G) ae eae 


Nearshore surface temperature gradient. 


2 


Mean coastal surface temperature. 


The maximum jet velocities that developed in experiments 1-4 exceed the 
observed values by about a factor of two. The exceptionally strong surface currents 
developed in experiments | and 3 may be explained in part by the unrealistically large 
value of wind stress at the coast in those experiments. The surface stress values were 
chosen by taking an offshore mean climatological value and applying it uniformly 
across the model domain. Experiments 2 and 4, which included the curl of wind stress, 
had more realistic (reduced) values of wind stress near the coast and indeed produced 
more realistic surface currents. Additionally, winds were input as steady which may 
produce stronger currents than a time varying forcing function. Time variable winds 
mav bv a Key to generating instabilities in the coastal jet and thus limiting its growth. 
Note that experiment 4 produced the maximum jet Velocities closest to those observed. 

Mince clish@remlocation of the jet axis, at the surface and with depth, compared 
favorably with observed locations. The vertical structure of the coastal jet tends to 
follow the coastal slope. The vertical eastern boundary in the model produced an 
almost vertically stacked jet, whereas the slope off northern California produces a 
coastal jet which tilts offshore with depth. Associated with the axis location is the 
offshore extent of the jet. Observations indicate that the coastal jet is a narrow feature 
which extends horizontally offshore from 25-100 km. Experimental results show a 
coastal jet about 100 km wide. The inclusion of wind stress curl (experiments 2 and 4) 
tended to increase the width of the modeled equatorward flow but the main body of 


the jet remained about 100 km in width. 


The depth of jet penetration in the water column is by far the largest deviation of 
miodel results from observations. Further experimentation and analysis is needed to 
identifv the cause of this major discrepancy. One possibilitv is the absence of a 
shelf’slope topography in the model. Model basin depth is constant at 4500 m which 
means that the coastal jet penetrated to about one-fourth of the basin depth. 
Observations show that the observed coastal jet penetrates to about one-third of the 
continental shelf depth. Therefore a relative comparison of model results with 
observations is good. Prelimunarvy results of a wind forcing experiment with idealized 
cross-shore topography indicates that the model coastal jet penetration is in fact 
dependent on shelf depth. An additional consideration in the vertical structure of the 
jet is the vertical temperature stratification. N{cCrearv er al. (1986) suggests that a 
shallow (~ 100 m) seasonal thermocline may be very important in the formation of 
the coastal jet, whereas the model’s basic stratification was that corresponding to the 
major thermocline having a vertical scale of 450 m. 

The poleward undercurrent maximum Velocities in experiments 1-4 were 
somewhat siow as compared with observations and also seem: independent of the 
magnitude of wind stress. Topography, again, may be an important mechanism in the 
formation of undercurrents. Open boundary conditions in the model may modify the 
formation of an undercurrent as compared with the effects of gyre scale circulation in 
the North Pacific. Wind relaxation events off the coast of California have also been 
considered to be closely associated with undercurrents and countercurrents within the 
CCS. Mooers ef al. (1976) stated that “The relationship between the northward 
Davidson Inshore Current, the California Countercurrent, and the undercurrent 
associated with coastal upwelling is not understood.” It is probably safe to say that 
that statement still holds true today. Recent emphasis in the oceanggragiie 
community has been placed on efforts to understand those relationships. The width of 
the undercurrent is generally considered very narrow (~ 15 km). as described bv Huyer 
and Kosro (1986). Experiment 4 produced the narrowest undercurrent (20 km) which 
corresponds well with observations. 

Vertical and horizontal shear measurements have demonstrated that the coastal 
jet off California and Oregon is a highly sheared system. Values for the vertical shear 
range from about 2 x 107s! to 5 x 103 s! (Huver, 1983: Kosro, 1986; Mooers e¢ 
ail., 1986). Experimental shear values fall within that range but on the low Sides 


ee o eaeine topographic influence may again be important to the establishment 
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of large vertical shear. The coastal shelf and slope, in compressing the vertical 
structure of the jet. may also increase the associated vertical shear. The measurements 
of horizontal shear in the CCS suggest that a tvpical value for the cyclonic shear may 
be on the order of f{,2 to f and that the anticyclonic shear is on the order of -f 5 
(Kosro, 1986). Experimental results demonstrate that a free-slip boundary condition 
On the eastern boundary (experiments | and 2) prevents the development of large 
horizontal cyclonic shear. Even with the zero-slip boundary condition (experiments 3 
and 4) the largest cyclonic shear developed by the PE model was .6f Anticvclonic 
shear was similarly smaller in the model experiments than the observed values of 
1) 3. 

Die Ose winewcOastal temperature structure is closely tied to the coastal jet 
structure owing to its baroclinic nature. The observed offshore surface temperature 
gradient during the upwelling season off northern California is on the order of .05 °:km 
with a mean Julv coastal temperature of 10-13 °C (Huvyer and Kosro, 1986). Model 
results show a range of gradients of .05-.07 ° km with coastal temperatures in the range 
of §.5-10 °C. Those values compare well, and the cooler model coastal temperatures 
mav be explained by the unrealistically high model jet velocities or inexact surface heat 
flux parameterization. 

It is difficult to make a quantitative comparison of experimental results with 
S@scrvalions Of the coastal jet in the CCS. Observed values of jet characteristics varv 
widely from study to study and the stable form of the classical two-dimensional coastal 
jet has been infrequently surveyed. From the above discussion, it is clear that the 
Structure of the simulated coastal jet is sensitive to the forcing function, boundary 
conditions, heating parameterization, topography, and a host of other conditions too 
numerous and complex for the scope of this study. Experiment 4, which included 
idealized, climatological wind stress curl with a zero-slip boundary condition on the 
eastern boundary, was the best simulation of the coastal jet of the four experiments 


conducted in this study. 
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VY. SUMIMARY AND CONCLUSIONS 


A. SUMMARY 

This study has examined the response of an idealized, flat-bottom, @asceme 
boundary oceanic regime to steady, equatorward, local wind-forcing on a seasonal time 
scale. Both uniform and offshore-varying wind, with zero-slip and free-slip boundarv 
conditions on the eastern boundary, were employed as forcing mechanisms in a ten- 
level, high resolution, primitive equation ocean model. 

A classical two-dimensional, upwelling-induced, equatorward flowing coastal jet, 
with an associated poleward undercurrent, developed in response to both umform and 
offshore-varying wind forcing. The variations in formation and structure of the coastal 
jet and undercurrent were explored and analyzed in four numerical experimieqme 
uniform wind with free-slip; offshore-varving wind with free-slip; uniform wind with 
zero-slip; and offshore-varving wind with zero-slip. 

1. Stability Analysis 

Instability of the mean flow did not occur in any of the aforementioned 
experiments and, therefore, eddies did not develop. A stability analvsis was conducted 
on all of the modeled coastal jets in this study along with an imposed coastal jet froma 
previous, unforced, flat-bottom experiment (Batteen er al., 1985) which formed eddies 
after approximately thirty days of model run time. Potential vorticity was computed 
throughout the model domain of each experiment and the cross-stream derivative of 
potential vorticity was found to change sign in the region of each modeled coastal jet. 
This fact established that the modeled coastal jets all satisfied the necessary condition 
for mixed (barotropic and baroclinic) instability. 

To further explore conditions for instability, a two-layer, quasigeostrophic. 
baroclinic instability model was applied to each experiment. It was established that all 
of the experiments satisfied the sufficient conditions for baroclinic instability as defined 
by the simple two-layer stability model. The four experiments conducted in this study 
did not significantly vary in critical wavelengths, most-unstable wavelengths, or in 
baroclinic instability e-folding times. However, the imposed baroclinic jet experiment 
(Batteen ef af, 1985) had a shorter critical wavelength and a shorter most-unstable 


wavelength than the four wind-forced experiments. The e-folding rates for longer 


78 


wavelengths were somewhat larger in the imposed barochnic jet experiment as 
compared with this study’s experiments; which would indicate that the wind-forced Jets 
are more unstable than the imposed jet at longer wavelengths. The reason whv 
instability of the mean flow did not occur in the model generated coastal jets remains 
an open scientific question. 

2. Comparison with Observations 

A comparison of model results with observations of the CCS was conducted. 
A numerical model can property simulate nature only to the extent which observations 
have provided the modeler with understanding of the feature he desires to simulate. 
This is a problematic area in modeling the CCS. Observations of the CCS demonstrate 
great variability in time and space and the classical two-dimensional upwelling-induced 
coastal jet has been infrequently observed and studied off California. With this in 
mind, a synthesis of literature addressing the coastal jet was conducted and a range of 
values for some of the more important characteristics was compiled and compared with 
model results. 

The model generated coastal jet compared favorably with observations in 
offshore axis location. undercurrent strength and width, vertical and horizontal shears 
and surface temperature gradients. Some deficiencies were noted. These included 
discrepancies in the strength, depth, and offshore extent of the modeled coastal jet 
compared with observations. Experiment 4, which included wind stress curl with a 
zero-slip boundary condition, 1s believed to be the most realistic model jet produced in 
this study. 

3. Effects of Offshore Varying Wind Forcing | 

Positive wind stress curl, within 200 km of the eastern boundary, was apphed 
by the offshore-varying wind forcing in two of the experiments (experiments 2 and 4). 
The positive wind stress curl induced cyclonic relative vorticity to the svstem and added 
an additional component of upwelling, due to the effect of Ekman pumping, to the 
coastal upwelling component (caused by equatorward wind forcing along the eastern 
niodel boundary). The Ekman pumping induced upward vertical motion, from surface 
divergence, within 200 km of the coast and created an extended zone of baroclinicity in 
the upper layers of the ocean. The extended zone of baroclinicity, therefore, increased 
the width of equatorward flow over that generated by the experiments with uniform 
wind forcing (experiments | and 3). Positive wind stress curl near the coast also moved 
the depth-averaged jet axis, along with the undercurrent axis, farther offshore than 


forcing without curt. 
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4. Effects of the Zero-Slip Boundary Condition 

As previously discussed, both free-slip and zero-slip boundary conditions were 
emploved on the eastern model boundary. The free-slip boundary condition implies 
that the offshore gradient of tangential velocity vanishes at the boundary (experiments 
| and 2), and the zero-slip condition sets the tangential velocity itself to zero at the 
coast (experiments 3 and 4). The zero-slip boundary condition, therefore, imparted 
cvclonic relative vorticity to the inshore side of the model coastal jet.” Prommemee 
comparison of model results with observations of the CCS, it is believed that eastern 
boundary coastal circulation ts best simulated, in a high resolution model, with the use 


of a zero-slip boundary condition on the eastern boundary. 


B. RECOMMENDATIONS 

This study was a first step in a numerical modeling effort at the Naval 
Postgraduate School in implementing wind forcing into a PE model designed to 
simulate the CCS off northern California. It is important for future studies to 
deternune the reasons why instabilitv of the mean flow produced by wind-forcing, in 
this study. does not occur when stabilitv analvses indicates that the flow is unstable. 
The development and evolution of eddies observed in the CCS is of extreme inteyese 
and understanding instability of the mean flow may be the key to understanding eddy 
processes. Designing future experiments to understand the vertical structure of the jet 
and why the model jet is much stronger and wider than the observed coastal jet should 
be given high priority as well. 

This study had the constraint of a flat-bottom: basin with a regular coastline. 
Implementation of topography first requires.the addition of depth-averaged flow within 
the PE model for scientific validity. Once the depth-averaged flow is incorporated, the 
door will be opened for a series of studies of wind-forced circulation over idealized 
topography. The influence of topography is thought to have a great impact on coastal 
circulation. Additionally, an irregular coastline which would simulate capes and bays 
along the coast of California would be a fascinating study. 

Finally, one should consider the inclusion of gvre-scale circulation as the baseline 
for wind-forcing experiments in the CCS. The gyre-scale circulation could be used to 
initialize the model prior to model spin-up bv local wind-forcing. Time-dependent 
winds could eventually be investigated as the prelude to model forcing by actual, 


observed wind records. 


SO 
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The time scale over which upwelling occurs may be analytically approximated 
through use of the principle of conservation of mass, scale analvsis, and Ekman 
transport. The upwelling rate can then be compared with the PE model results for 
verification of model physics. The approximation involves no diabatic heating and 
Meeglects the Component of velocity due to pressure gradient force. A control 
Experiment was conducted which included no diabatic processes (surface thermal 
damping) and which had uniform, steady wind forcing. That experiment will be used 
for comparison with the following approximation. All svmbols. in the following 
derivation, are as defined in Chapter II unless otherwise stated. 

Momarmive 20 a scale for the vertical comipoment of velocity, the equation of 


continuity, with the assumption of incompressibility, 1s written as 


Gu OV Ow 
Ox Ov Oz 





=O (A.1) 


Assuming no longshore variation in v, then equation (A.1) may be reduced to 


Gu Ow | 
a= - == Pero . (A.2) 
Ox OZ 


Scale analysis of equation (A.2) vields the following scale for w (vertical velocity): 
Cae 
E 


X 


Wh es 





; (A.3) 


where UH is mass transport and L, is a length scale. 
Since it is desired to isolate the effect of upwelling caused only by motion of the 
fluid away from the solid boundary driven by surface stress, Ekman miass transport 


(Mp) is required where 





Mise = (eae pe Up dz = ; 


which may be rewritten as 
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OF 


LH = : , (ee) 





Using equation (A.3), W 1s substituted into equation (A.4) to arrive at the 
following relationship between vertical velocity and surface stress: 
T 


W ~ ~ ; (A.5) 
D, ees 





The surface stress can be calculated, through use of the bulk aerodynanuc formula with 


a constant drag coefficient, by 


tT = ~p,Cp ¥° 


Vv 


Substituting this expression for e into equation (A.5) results in the final form for 


vertical velocity (W) of 


Pp, op 
p, FL 


wf (4.6) 


xX 
Using a value of V = 900 cms‘! and L. = 14.2 km (the second internal Rossby radius 
of deformation computed by the method of Feliks (1985)). equation (A.6) vields a scale 


value for the vertical velocitv, generated by coastal upwelling, of 
W =~ 1.00947 cm, see. =" S715 muda. 


Model output demonstrates that. in the upper levels of the ocean, vertical 
velocity decreases by 80-90% within two gridpoints of the coast. This length scale, of 
16 km, compares Well with the second internal Rossby radius of deformation. The 
model length scale of vertical velocity does not, however, compare well with the first 
Rossby radius of deformation which was computed to be 29.3 km. The horizontal 
length scale corresponding to the second Rossby radius of deformation 1s appropriate 
because the model vertical velocity has a maximum close to the sea-surface. 

A comparison of the rate of model upwelling with that obtained from an Ekman 
balance was also made starting. arbitrarily. at model day 30 since Ekman balance ts not 
well established during the initial days of model spin-up. At this time, the depth of the 


9.03° C isotherm (chosen for illustrative purposes) at the coast was 169.4 m (level 4). 


Using the vertical velocity approximation computed above, the 9.03° C isotherm 
Simi aeee displaced to level 1 (12.9 m) in 20.65 days. \lodel output demonstrates that. 
pea 50-G))icOastal temperature at level lis 9.19" C which compares quite favorablv 


with the analytically approximated temperature of 9.03° C. 
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